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ABSTRACT 

Geometrical cross sections of dust aggregates determine their coupling with disk 
gas, which governs their motions in protoplanetary disks. Collisional outcomes also 
depend on geometrical cross sections of initial aggregates. In the previous paper, we 
performed three-dimensional iV-body simulations of sequential collisions of aggregates 
composed of a number of sub-micron-sized icy particles and examined radii of gyration 
(and bulk densities) of the obtained aggregates. We showed that collisional compression 
of aggregates is not efficient and that aggregates remain fluffy. In the present study, we 
examine geometrical cross sections of the aggregates. Their cross sections decreases due 
to the compression as well as their gyration radii. It is found that a relation between 
the cross section and the gyration radius proposed by Okuzumi et al. is valid for the 
compressed aggregates. We also refine the compression model proposed in our previous 
paper. The refined model enables us to calculate the evolution of both gyration radii 
and cross sections of growing aggregates and reproduces well our numerical results of 
sequential aggregate collisions. The refined model can describe non-equal-mass collisions 
as well as equal-mass case. Although we do not take into account oblique collisions in 
the present study, oblique collisions would further hinder compression of aggregates. 

Subject headings: ISM: dust, extinction — planetary systems: protoplanetary disks 



1. INTRODUCTION 

In the theory of planet formation, planets are thought to have formed in protoplanetary disks 
through mutual collisions and coalescence of planetesimals. The formation process of planetesimals, 
on the other hand, still has a large uncertainty. Before the planetesimal formation, dust grains grow 
through their collisional coalescence in a protoplanetary disk and settle to the disk mid-plane, which 
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form s a dense dust layer at the mid -plane (e.g.. ISafronovlll969l : iNakagawa et al.lll98ll : iTanaka et al 
20051 : iDullemond &: Dominik |2005|). Planetesimals would be formed in the dust lay er through 
gravitational ins tability (e.g., Goldreich &: Ward 1973 : Sekiya 1998^:Vo udin Shull2002 ). streaming 



instability (e.g.. lYoudin &: Goodman 



Weidenschilling Cuzzil 119931 : 



2005 



Brauer et al 



Youdin fc Johansen .2007i ) or simple coalescence (e.g., 



2008a 



In these models of planetesimal formation, 
motion of dust grains is an important factor because it determines the spatial distribution of dust 
grains and the collision speed between them. Furthermore, their motion is governed by the drag 
forces from the disk gas. 

Gas drag forces on dust grains strongly depend on their internal structure (or their bulk den- 
sities). Most studies of dust growth in protoplanetary disks have assumed compact structure of 
dust grains. However, dust grains growing through mutual collisions would actually be aggregates 
of (sub-micron) primitive grains and the aggregates have a fluffy structu re with an extremely low 
bulk density, as reported by experimental and theoretical studies (e.g. , iBlum k. Schrapled |200J; 



Ormel et al.l 120071 : ISuyama et al 



2008 



Okuzumi et al. 



2009 



Zsom et a] 



2011 



Such fluffy aggre- 
gates have large ratios of their geometrical cross sections to masses, which significantly enhance gas 
drag forces on them compared with compact dust grains. Hence, in order to clarify dust growth 
and planetesimal formation in protoplanetary disks, we have to examine the internal structure 
and the geometrical cross sections of dust aggregates (we use a term 'cross section' in referring to 
'geometrical cross section', hereafter). 



Suyama et al. (j2008l ) (hereafter S08) performed A^-body numerical simulations of sequential 



aggregate collisions to examine the compression process of growing aggregates. The sequential 
collisions mean that we repeat collisions of aggregates obtained at the previous collisions. With such 
a simulation, we can observe a natural evolution of the aggregate structure. Their numerical results 
showed that large aggregates have an extremely low bulk density in spite of compression at aggregate 
collisions. In the early stage of dust growth, aggregates just stick without any restructuring because 
of their low impact energy and they have a fluffy structure with an extremely low bulk density as 
they grow. In the later stage in which the impact energy exceeds a critical energy, aggregates are 
gradually compressed. Even in this compression stage, their density remains very low. It is found 
that the compressed aggregates have a low fractal dimension of 2.5. This structural feature causes 
the low density of the compressed aggregates. SOS also derived a formula describing the density 
evolution of growing aggregates. To estimate their bulk densities, S08 used the so-called gyration 
radii of the aggregates but did not examine their cross sections. However, cross sections are directly 
related to the gas drag forces rather than gyration radii. It is necessary to clarify the evolution of 
cross sections of dust aggregates during their growth. 

Cross sections of aggregates depend on their internal structure. There are two simple aggre- 
gate models. One is Ballistic Cluster-Cluster Aggregation (BCCA). A BCCA cluster is formed 
through collisions between two equal-sized clusters. Second is Ballistic Particle-Cluster Aggrega- 
tion (BPCA). A BPCA cluster is formed through deposition of small monomer particles on a large 
cluster. For both BCCA and BPCA, restructuring is assumed to be negligible at each collision. The 
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BCCA clusters have very fluffy and open structures and the BPCA clusters have relatively compact 
structures. Figure [1] shows the ratios of cross sections to masses of aggregates. It is shown that 
the cross section per mass strongly depend on aggregate types. Cross sections of dust aggregates 
are expected to be between those of BCCA and BPCA clusters. One may consider that a cross 
section is approximately given by the square of a gyration radius. Cross sections of aggregates are, 
however, generally independent of their gyration radii, especially for highly fluffy aggregates. The 
non-dimensional ratio of the cross section to the square of the gyration radius gradually decreases 
with their growth for BCCA clusters (Minato et al. 2006; see also Fig. 6). On the other hand, this 
ratio is almost constant in the growth of BPCA clusters. Okuzumi et al. (2009) proposed a useful 
relation between the cross section and the gyration radius for various aggregates formed through 
hit-and-stick growth (as well as BCCA and BPCA). Paszun and Dominik (2009) also derived an- 
other relation. Nevertheless, it is not clear whether these relations are also valid for aggregates 
compressed at collisions. We check the validity of these relations, using the resultant aggregates 
obtained by SOS. 

Once we find a valid relation between the cross sections and the gyration radii, it would be 
very helpful to describe the evolution of the cross sections because the compression model by SOS 
can describe gyration radii of growing aggregates. The compression model by SOS, however, has 
some limitations. This model is not directly applicable to low-energy collisions (i.e., hit-and-stick 
collisions) or to non-equal-mass collisions. In order to describe gyration radii and cross sections of 
aggregates for all growth stages seamlessly, we further refine the compression model, by removing 
these limitations. 



In laboratory experiments, IWeidling et al.l (j2009l ) examined compression of aggregates consist 



ing of 1.5/im-diameter Si02 spheres at their multiple rebounds and also developed an empirical 
compression model. Initial aggregates in their experiments possess a volume filling factor of ~ 0.1, 
which is approximately equal to that of BPCA clusters. On the other hand, SOS and the present 
paper focus on the compression of fluffier aggregates of which filling factor is between BCCA and 
BPCA clusters during their collisional growth. Hence our compression model and theirs are com- 
plementary to each other. As mentioned above, dust aggregates are expected to have much smaller 
bulk densities than BPCA clusters at the early stage of their growth in protoplanetary disks. Our 
compression model is useful as long as bulk densities of aggregates is lower than that of BPCA 
clusters. 

It should be noticed that SOS only considered head-on collisions of aggregates in their numerical 
simulations of sequential collisions. The oblique collisions are expected to hinder the compression 
(Wada et al. 2007; Paszun and Dominik 2009). In the present study, however, we use the results 
in SOS as the first step. We will examine the effects of oblique collisions in future work. 

In the next section, we briefly summarize the results of SOS. In Section 3, we numerically 
calculate cross sections for the aggregates obtained by SOS. We find that Okuzumi et al's relation 
between the cross section and the gyration radius is valid for compressed aggregates, too. In Section 
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4, we refine the compression model by SOS, by removing its limitations in a reasonable way. We 
find that the refined compression model reproduces well both of gyration radii and cross sections 
of aggregates obtained by the numerical simulation, with the help of Okuzumi et al's relation. We 
also check the validity of the refined model for non-cqual-mass collisions with additional numerical 
simulations of aggregate collisions. A summary is given in the last section. 



2. RESULTS OF AGGREGATE COMPRESSION IN AT-BODY SIMULATIONS 

BY SOS 

Suyama et al. (2008) performed A?^-body numerical simulations of head-on aggregate collisions 
and examined the density evolution of aggregates growing through the collisions. Wc examine the 
cross section of the resultant aggregates obtained by SOS. Before that, we briefly describe the 
numerical results of SOS. 

In the simulations, aggregates consist of a large number of icy spherical particles with a radius 
of ri = O.l/im. SOS adopted the particle-interaction model by Wada et al. (2007). In the interaction 
model, repulsive and adhesive forces in the normal direction between particles in contact are given 
by the JKR theory (Johnson et al. 1971). A tangential force and a torque also arise to resist the 
slide, roll, and twist motions between them. Aggregate compression is regulated mainly by inelastic 
rolling motions of the constituent particles (e.g., Dominik &; Tielens 1997; Wada et al. 2007, 200S 
[hereafter W07,W08] ; Giittler et al. 2010). The rolling energy -Broil is the energy required for 
rolling a particle on its contact neighbor by an angle of 7r/2. The rolling energy is given by (W07, 
SOS) 

Eron = 67r^7ri^crit, (1) 

where ^crit is the critical displacement for inelastic rolling motion. The parameter range of ^crit is set 
to be from 2 to 16 A in SOS. A large rolling energy E'rou suppresses the restructuring of aggregates. 
To examine the structure evolution of growing aggregates, SOS performed AT-body simulations 
of sequential collisions. Each simulation starts from a collision of aggregates composed of two 
particles (i.e., dimers) and ends with a collision of aggregates composed of 16,384 particles. The 
resultant aggregate obtained in the previous collision is used as initial aggregates at each collision 
in the simulation of sequential collisions. The impact velocity is constant in sequential collisions. 
For various (constant) impact velocities and critical rolling displacements, they performed a large 
number of runs of sequential collisions. 

As an index of the size of an aggregate, SOS adopted the radius of gyration, rg, defined by 



N I 19 



iV ' 



where Xi is the position of particle i, xm is the position of the center of mass of the aggregate, and 
N is the number of particles composing the aggregate. Using the radius of gyration, the volume V 
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and the bulk density p of the aggregate are evaluated to be (Mukai et al. 1992; W08) 

respectively, where y^S/ Sr^ is the so-called characteristic radius of an aggregate and mi is the mass 
of a constituent particle. 

Aggregates are expected to have a BCCA structure for collisions at sufficiently low velocity 
because of sticking together of equal-mass aggregates without any restructuring. If the compression 
is effective at collisions, the gyration radii of aggregates would become smaller than those of BCCA 
clusters. It is meaningful to compare the obtained aggregates with BCCA clusters. Since BCCA 
clusters have a fractal dimension of ~ 2, the radius of gyration of the BCCA cluster is given for 
large N by (e.g., Mukai et al. 1992; W08 )Q 



r-5,BCCA ^ iV°-''°ri. (5) 

Figure [2^ shows the gyration radius of the aggregates in the simulations of sequential collisions 
performed by SOS for various values of parameters, ^crit and the impact velocity fimp- The density 
of monomer particles is given by Pm{= 3mi/[47rri]). In the simulation of sequential collisions, the 
size of growing aggregates is dependent on the direction of each collision. SOS did 30 runs of the 
simulation of sequential collisions and obtained the averaged value of rg from 30 runs for each ^crit 
and Uimp- In Figure [2^, the horizontal axis is the number of the constituent particles, A'^, in the 
growing aggregates and the vertical axis is the gyration radii divided by N^/'^ri for comparison with 
BCCA clusters. The dashed line represents the radius of the BCCA cluster and it is almost flat for 
large as expected from equation ([5]). The size of small aggregates produced in our simulation 
is almost the same as that of BCCA clusters. This is because the impact energy is small enough 
at the early stage of the aggregate growth and the compression is ineffective at each collision. As 
the aggregates grow, the impact energy increases. When the impact energy attains to -Broil, the 
compression of the aggregate starts: aggregates become smaller than BCCA clusters. The critical 
number of particles, A'crit, in the aggregate for compression is given by (WOS, SOS) 

A^cnt = /3^%^, (6) 



^Exactly speaking, equation (O is satisfied for BCCA clusters formed through head-on (hit-and-stick) collisions, 
which have the fractal dimension of 2.0. When offset collisions are also included at the formation of BCCA, their 
fractal dimension is 1.9 and the gyration radii proportional to A^"'^^ (Okuzumi et al. 2009). Although the later 
BCCA is more realistic, the former BCCA is used in SOS and the present study since SOS consider only head-on 
collisions in their simulations. 
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where /3 is a non-dimensional coefficient. In Figure [2^, we also plot the critical number A'crit with 
filled circles on each curve, by setting /3 = 0.5. Figure[2]3 shows bulk densities of growing aggregates 
in these simulations, which evaluated with Vg by equation @. We also plot the density of BCCA 
clusters. It decreases as N~^'^^ for large A^. After the onset of compression (i.e., N > Acrit); the 
bulk densities of the aggregates obtained by the numerical simulation are larger than that of the 
BCCA but still keep on decreasing gradually in all cases, indicative of the inefficiency of collisional 
compression. SOS also developed the compression model, which reproduces the density evolution of 
growing aggregates in the compression stage. We will describe the compression model in Section 4. 



3. GEOMETRICAL CROSS SECTIONS OF AGGREGATES PRODUCED IN 

THE SEQUENTIAL COLLISIONS 

We numerically calculate cross sections of resultant aggregates obtained by SOS. The cross 
section of a dust aggregate is given by the area of the shadow of the aggregate projected onto 
a plane. The area of the shadow is calculated by counting the number of square meshes in the 
shadow (Fig. [3]). The width of the square meshes is set to be 0.0055ri. This width is much smaller 
than the radius of the monomer particle, ri, though meshes with a much wider width are drawn 
to emphasize them in Figure [3l The area of the shadow is dependent on the plane onto which the 
shadow is projected. We calculate the areas of the shadows for 30 orientations randomly chosen and 
define the cross section of the aggregate by the mean values of the areas. Figure U] shows the cross 
section calculated in this way for aggregates produced in the simulations of sequential collisions. 
The vertical axis is the cross section divided by Nnrf, which corresponds to the (non-dimensional) 
cross section per mass. If the overlapping of the monomer particles in the shadow is negligibly small, 
the value of the vertical axis approaches unity. A filled circle in Figure \^ indicates the shadow area 
for each orientation and the line shows the mean of them. The cross section per mass decreases as an 
aggregate grows in the simulation of sequential collisions due to the overlapping of the constituent 
particles. Since SOS did 30 independent runs of sequential collisions, the mean cross sections are 
calculated and plotted as thin lines in Figure Hb. Then we obtain the averaged value of the 30 
mean cross sections as shown by the thick line in Figure Hb. The dispersion of the mean geometrical 
cross section can be evaluated in Figure SJ^. The standard deviation of the mean geometrical cross 
section is equal to or less than 11% of its averaged value during the aggregate growth. In this way, 
we did two kinds of averaging to calculate the cross section of growing aggregates in the simulation. 
We show and discuss the cross sections of aggregates by using these finally-obtained cross sections 
hereafter. 

Figure [5] shows the cross section of resultant aggregates for various values of the parameter set 
(^criti ^imp)- The averaged cross section of the BCCA clust er is also calculated and plotted. The 



cross section of BCCA we obtained agrees with the result of lMinato et al.l (j2006l ). Even for BCCA, 
the ratio S/{NTTr1) gradually decreases with an increase in N due to the overlapping of constituent 
particles. In the early stage of the aggregate growth (i.e. for small N), the cross sections of the 
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resultant aggregates change almost along the line of BCCA. As the aggregates grow, however, their 
cross sections deviate from the line of BCCA and become much smaller than that of BCCA, which 
is due to compression at collisions. This qualitative tendency is consistent with the evolution of the 
radius of gyration (Fig. [2^). Although the change in the cross sections is gradual and the starting 
points of compression in the cross sections are not clear compared with those in the gyration radii, 
the starting points are also described by equation ([6]), by setting the parameter /3 to be 2.0. The 
larger value of (3 than Figure [2] indicates that the onset of compression in the cross section is later 
than that in the gyration radius. 

In Figure [6] we plot the ratio of the cross section S to vrr^ for all resultant aggregates with solid 
lines. The ratio of S to irr'^ decreases for small aggregates, which is consistent with the BCCA case 
(dotted lines). For sufficiently large aggregates, the ratio increases as a result of their compression. 
Okuzumi et al. (2009) proposed a useful expression of the cross section S for aggregates formed 
through hit-and-sticks. The expression is given by 



Equation ([5]) is used as the expression of rg,BCCA(A^)- We also plotted the cross sections obtained 
from equation ([7]) with dashed lines in Figure El It is found that the expression by Okuzumi 
et al. (2009) reproduces the cross section surprisingly well for compressed aggregates as well as 
hit-and-stick aggregates, by using the gyration radius r^. In this expression, the information of 
compression is correctly included through the gyration radius. Paszun and Dominik (2009) also 
derived another relation between S and the aggregate size (i.e., eq.[ll] of their paper). Figure 7 is 
the same as the left-bottom panel of Figure 6 but the prediction by Paszun and Dominik is also 
plotted. Although the prediction by Paszun and Dominik is consistent with the numerical results, it 
overestimates S when the ratio S/{'Kr'^) is larger than unity (i.e., for relatively compact aggregates) 
and underestimates for 5/(7rr^) < 0.7. The underestimation was also reported by Okuzumi et al. 
They found that the underestimation in Paszun and Dominik's model is severe especially for large 
and fluffy aggregates. For other we also find the same trend as in the case of = SA shown 
in Figure [3 Hence it is concluded that the model of Okuzumi et al. is more accurate than that of 
Paszun and Dominik. Once an accurate compression model describing Vg is obtained, it enables us 
to calculate the evolution of the cross section with the help of Okuzumi et al's model. 



1 




(7) 




(8) 



4. 



COMPRESSION MODEL 



A compression model describing rg was developed by SOS but it has two limitations. The 
model is not directly applicable to low-energy collisions (i.e., hit-and-stick collisions at the early 



- 8 - 



growth stage) or to non-equal-mass collisions. In order to describe gyration radii and cross sections 
of aggregates for both the early hit-and-stick stage and the compression stage seamlessly, we refine 
our compression model, by removing these limitations in a natural way. Before that, we briefly 
describe the compression model by SOB. 



4.1. Compression Models of W08 and SOS 

W08 developed a compression model by introducing the pressure (or the strength) of aggregates 
to explain their numerical results on collisions between BCCA clusters. The compression model of 
SOS is based on that of W08. At a collision of two aggregates with the impact energy -Bimp, the 
compression of the merged aggregate from the initial volume, Vinitiah to the final volume, Vfinah is 
described in the model of W08 by 

/■Vflnal 

^imp = - / PdV. (9) 

^initial 

The initial volume initial is defined by the volume of the merged aggregate at the moment that the 
two aggregates just stick. After the moment of the sticking, the compression proceeds. The volumes 
before and after the compression are evaluated with the radius of gyration, Vg, as in equation 
The pressure P of the aggregates is given by 

P = 2 ( ^-Y ^-^^ ( (10) 



3y mi \Pm) 

where the fitting parameter b is set to be 0.15. Note that the pressure P of the aggregates is 
dependent on the total number of constituent particles (or the total mass) as well as the density. 
That is, P is not an intensive variable. This strange property in the pressure comes from the 
fractal structure of the aggregates. W08 showed with their simulation of collisions between BCCA 
clusters that the compressed aggregates have internal structures with a fractal dimension of 2.5. 
The simulation of sequential collisions done by S08 showed that their resultant aggregates also have 
the same fractal dimension of 2.5. 

In order to describe the compression of such fractal aggregates, W08 also introduced the fractal 
volume defined by 

Vfirg) ^ arl\ (11) 

where the coefficient a is given by (97r/5)^/'^r(9/4) ~ 7.7. Using the fractal volume, the fractal 
density is defined by 

„ir,) . ^ = ^rf \ (12) 



The dimensions of the fractal volume and the fractal density differ from those of the ordinary 
volume and density. These fractal quantities are related with the ordinary quantities V and p as 
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3/2 r^f^ N\6/5 



] f N-y\ (14) 



Pm \5J V Pf,l 

where Vm (= 4:/37rrf) is the volume of a monomer and Vf^i is given by arf'^. Using the fractal 
volume Vf (and the fractal density pj), equation @ is rewritten as 

/■V/, final 

Eimp = - PfdVf, (15) 

■'Kf, initial 

where the fractal pressure Pf is given by 

P^^P^ = ^^Ji^(Il)\ (16) 
dVf mi Vp/,1/ 

It should be noticed that the fractal pressure is dependent on pj but not on the total mass. That 
is, Pf is an intensive variable. Equation (jl5p (or [9]) reproduces the numerical results on the 
compression at collisions between BCCAs. 

SOS pointed out that WOS's compression model needs a minor modification to describe the 
compression of partially compressed aggregates at their collisions, which occurs in their simulation 
of sequential collisions. At the moment of sticking at each collision, large voids are produced in 
the merged aggregate. The volume of the new voids is included in the initial volume of the merged 
aggregate, Vj^initiai in equation ([15]). The energy required for compression of the new voids is 
~ -Eroii and it is much smaller that that predicted by equation (jl5p at collisions between partially 
compressed aggregates. To describe the compression at such collisions, SOS modified WOS's model. 
Since the energy required for the crush of the new voids is negligible, the initial fractal volume of 
the merged aggregate in equation ([TS]) is set to be the sum of the fractal volumes of two colliding 
aggregates, by removing the volume of the new voids. That is, 

V/,initial = +y/,2. (17) 

Using equations (fT5]) - p!7|) . we have the (final) fractal density of the merged aggregate, Pffinaiy 
produced at collisions of two equal-mass aggregates with the fractal density pj^ (= Nmi/[Vf^i + 

\ 4 / \ 4 CI 
P/,final \ _ I Pf,0 \ ^ -^imp ^^^-^ 



Pf,l ) \Pf,l) bNEj-oU 

where N is the number of constituent particles in the merged one. Equation (jlSh describes the den- 
sity evolution of partially compressed aggregates growing through mutual collisions. Equation (jlSp 
with b = 0.15 reproduces the density evolution of growing aggregates in Figure [2] for > Ncrit, as 
seen in Figure S of SOS. 



4.2. Refinement of the Compression Model 



The compression model by SOS is not applicable to the early growth stage {N < Acrit) where 
the compression is ineffective. This is because the energy required for the crush of the new voids is 
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neglected in the model. This limitation is removed, by taking into account the compression process 
of the new voids produced at sticking of two aggregates. Furthermore, the model of SOS assumes 
the collisions of equal- mass aggregates. We also remove this limitation in a reasonable way. At the 
extension of our compression model to non-equal-mass collisions, we assume that the compressed 
aggregates have the fractal dimensions of 2.5 as well as in the case of equal-mass collisions. The 
validity of the assumption will be discussed in Section 5. 

Density evolution of aggregates at each collisions is divided into the following three steps: 

1. Creation of new voids at the sticking of two aggregates. 

2. Compression of the new voids in the merged aggregate. 

3. Further compression of the merged aggregate after the crush of the new voids. 

These steps are schematically explained by Figure El We describe the change in the volume or the 
density of the aggregates at each step in detail below. 

In Step 1, the density decreases because of the new voids in the merged aggregate. The change 
in the density is described by Okuzumi et al. (2009). When the two aggregates with the masses Mi, 
M2, and the volumes Vi, V2 (< Vi) collide with each other, the volumes of the merged aggregate 
^1+2 (before the compression) is given by 

K+2 = Vi+V2 + Koid, (19) 

where the volume of the voids Vyoid is obtained from the empirical formula 



Koid = min 



XBCCA-1.031n( ^^^^^ ] ,6.94 



V2 (20) 



(see Okuzumi et al. [2009] for the derivation). The collision of equal-mass aggregates is not assumed 
in Okuzumi et al. (2009). Note that xbcca(= 2^/"^^ - 2) is 0.83 for BCCA clusters formed through 
head-on collisions having the fractal dimension dj = 2.0. The density of the merged aggregates 
before the compression is given by (Mi + M2) /V(_^2- "^^^ fractal volume ^1+2 fractal 
density 1+2 obtained from equations (fT3]) and ([HI), respectively. 

In Step 2 and 3, the merged aggregate is compressed and the density increases. Step 2 is the 
compression of the new voids. We put the energy required for the crush of the new voids to be 
6'-Eroii, where b' is the non-dimensional parameter. By fitting with results of the simulation, this 
parameter is fixed to be b' = 3b, as will be shown in the next subsection. This relation would be 
reasonable because both parameters are related to the beginning of compression. 

When the impact energy -Eimp is smaller than b'Ej-oii, the new voids are only partially com- 
pressed at Step 2. In this case, the compression at Step 3 does not occur because of the small 
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impact energy. The impact energy -Eimp is evaluated using the reduced mass, (1/Mi + I/M2) ^. In 
this case, the decrease in the fractal volume of the merged aggregates, AVj i_|_2; is evaluated to be 

= T^^^/.void, (21) 

where 

^/,void = 17',l+2-^/,l-^/,2. (22) 

It should be noticed that the fractal volume Vj^void defined by equation (j22p is negative in collisions 
where the volume ratio V1/V2 is larger than 8 x 10'*. A prescription for such high- volume-ratio 
collisions will be described at the end of this subsection. Here we consider the case where Vj^void is 
positive. Then the final fractal volume Vf^i^2 after the compression is given by 

Vf,l+2 = y'f, 1+2 - A^/,1+2. (23) 

and the final fractal density /9/,finai of the merged aggregate is given by (Mi + M2) /Vf,i+2- 

When the impact energy is larger than b'Ej-on, the new voids are crushed completely and the 
merged aggregate is further compressed at Step 3. Since the fractal volume of the merged aggregate 
is Vf^i + Vf^2 at the end of Step 2, the final fractal volume after Step 3 is obtained from the 

equation 

^imp - b'E,„u = - f Pf{Pf)dVf. (24) 

Integrating the RHS of equation we obtain the final fractal density for i?imp > ^'-E-roU as 

/'/.finaA^ _ / P/,o\^ -E'imp — b' Ej-^n 



+ — ^ , (25) 

where we used 

M1 + M2 

In the limit of -Bimp ^ ^'-E'roU) equation (p5]) is identical to equation (fTSj) (or the compression model 
by SOS). The final fractal volume Vy,i+2 is given by (Mi + M2)/ Pf final- Z,From the fractal density, 
we obtain the gyration radius Vg, using equation ()12p . The cross section S is also obtained from 
Okuzumi et al's expression (eq.[7])- In this way, we can calculate the density evolution (i.e., the 
evolution of rg and S) at both low- and high-energy collisions, using equation (p3|) for -Eimp < b'E^all 
and equation (|25]) for E'imp > ft'^roii- 

In high- volume-ratio collisions where V1/V2 > 8 x 10*, as noticed above, Vy,void is negative 
and a special prescription is necessary. Since a negative V/^void means no voids. Step 2 should be 
omitted and the merged aggregate is compressed only with Step 3. That is, equations (pl |) - ([26|) 
are used for all impact energies in this case. In equations ()24p -()26p. the terms of b'E^oW is omitted 
and Vf^i + Vf^2 is replaced by because Step 2 does not occur. 
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4.3. Test of the Refined Compression Model 



Let us test the refined compression model with the numerical results. Using the refined com- 
pression model, we calculate the evolution of gyration radius for the same condition as the numerical 
simulations by SOS and also obtain the cross sections of the aggregates with equation ([7|). The re- 
sults are shown in Figures [9] and [101 In Figure [U we plot the evolution of the gyration radius 
calculated with the refined compression model and compared it with the numerical results by SOS. 
The parameters b and b' are set to be 6 = 0.15 and b' = 36 (= 0.45), respectively. With this setting 
of the parameters, the refined model reproduces well the numerical results at both the early growth 
stage and the compression stage. Figure [TOl shows the evolution of the cross sections and indicates 
that the refined model also succeeds in describing the cross sections with the help of equation ([7]). 

In the above, the evolution of gyration radius of growing aggregates are calculated with the 
refined compression model and their cross sections are indirectly calculated, by using rg and equa- 
tion ([7|). We also propose another way to describe the cross sections of aggregates. We define 
alternative characteristic sizes of aggregates rs by 



It would be possible to describe the evolution of rs directly (instead of the gyration radius) with 
the refined compression model in the following way. Using this characteristic size rs instead of 



Y^5/3rg, we can define the volume and the bulk density of the aggregate by the similar equations 

to ([3]) and (j4j). The fractal volume and the fractal density are also defined in the same way. Then, 
applying the refined compression model to the fractal density defined with rs, we can describe the 
evolution of rs as well as in the case of rg. The evolution of the cross section S is calculated with 
equation ()27p . This is a direct way to describe the cross section rather than the above. In this 
calculation of S, we have to be cautious with the following two points. One is the modification 
in equation ([T9|). At a collision of sufficiently fluffy aggregates, equation (fTO]) can give the volume 
U/_,_2 larger than that of the BCCA cluster with the same mass, T^ccA when the size rs is used 
instead of rg. Such a large V(^2 is not realistic. In this case, we set the volume V(_^_2 = Vbcca 
instead of equation (|19p . The other point is the parameter b. Although b is set to be 0.15 in the 
case of rg, we have to calibrate the parameter b again in the case of as a result of the fitting 
with numerical results. In Figure \TT\ we plot the evolution of S calculated with this direct way. In 
this calculation, the parameters are set to be 6 = 0.6 and b' = 3b. We see that the refined model 
works well for the evolution of the cross sections with this direct way, tool 

It is found that the refined compression model enables us to describe the whole evolution of 
the radius of gyration rg and the cross section S of growing aggregates. Note that the refined 
model is applicable to non-equal-mass collisions though the tests in Figure [9lfTT] are done only 



rs = y^SjiT. 



(27) 




In the early growth stage of Figure 1111 we used V{^2 = Vbcca instead of equation (|19l) when the volume of 
equation (|19p is larger than Vbcca, as mentioned above. 
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in the equal-mass case. At the extension of the refined model to non-equal-mass collisions, the 
fractal dimension of compressed aggregates is assumed to be 2.5 even in the case of non-equal-mass 
collisions though it is not verified with A^-body simulations in non-equal- mass cases. It is possible 
that a very large mass ratio increases the fractal dimension, as seen in BPCA clusters. Okuzumi 
et al. (2009) examined the effect of the mass ratio on the fractal dimension for aggregates growing 
with hit-and-stick collisions with A^-body simulations. They showed that the collisions with the 
mass ratio of 10 increases the fractal dimension df only by 0.1 (see their figure 6). Furthermore, 
collisions with such a mass ratio have a major contribution at dust growth in protoplanetary disks, 
as shown by Okuzumi et al. (2009,2012). Hence the effect of non-equal-mass collisions would not 
change largely the fractal dimension of compressed aggregates in the realistic growth process. 

To confirm the validity of the refined model in non-equal-mass collisions, we have further 
performed additional A^-body simulations of aggregate collisions. Similar to W08, we consider 
collisions of two BCCA clusters but their masses are not equal in the present case. The projectile 
BCCA cluster consists of 1024 particles (or 4096 particles) while the number of constituent particles 
of the target BCCA cluster is 16384. Their mass ratio is 1/16 (or 1/4). The constituent particles 
are icy ones with the radius with 0.1/i m. The impact velocity fimp is a parameter. We set 
■Wimp < 4.4ms~^ since we focus on the compression process rather than fragmentation (Wada et 
al. 2007,2008). The numerical results of compression at the non-equal-mass collisions are shown 
in Figure [T2l The predictions by the refined model with b = 0.15 are plotted by solid lines and 
the dashed line indicates the formula by W08 (their equation [45]). The numerical results in the 
non-equal-mass collisions approximately agree with the predictions by the refined model though 
upper shifts by ~20% are observed in the case of Mp/Mt = 1/16. As well as in the case of equal- 
mass collisions, the slope in rg-E^^p relation is approximately given by -0.1, which indicates that 
compressed aggregates have the fractal dimension of 2.5. The upper shifts in the numerical results 
indicate that the compression requires larger impact energy in non-equal-mass collisions than the 
equal-mass case. This effect in non-equal-mass collisions would be included by adopting a larger 
parameter b in the refined model. Furthermore, for high-mass-ratio collisions with Mt/Mp ^ 10, 
the refined model is not verified through A^-body simulations yet although such collisions have only 
a minor contribution in dust growth (Okuzumi et al. 2009,2012). In the future work, the effect 
of non-equal-mass collisions should be further examined in the numerical simulation of sequential 
collisions as done by SOB in order to calibrate the parameter b more accurately. 



5. SUMMARY 

We examined the evolution of the geometrical cross section of the growing (icy) aggregates ob- 
tained by A^-body simulations of sequential head-on collisions (S08) and constructed to construct 
a refined compression model, which is applicable to the description of the evolution of both geo- 
metrical cross sections and gyration radii of growing aggregates. The results are summarized as 
follows: 
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1. We examined geometrical cross sections of the aggregates produced in the simulation of se- 
quential collisions done in our previous paper. As aggregates grow, compression becomes 
effective and makes their cross sections smaller than those of the BCCA clusters. The be- 
ginning of the compression is given by equation ([6|), as seen in the evolution of the gyration 
radius. 

2. The relation between the cross section and the gyration radius seen in aggregates obtained 
by SOS is well described by Okuzumi et al's expression. This indicates that Okuzumi et 
al's expression is valid for compressed aggregates as well as hit-and-stick aggregates. If the 
evolution of the gyration radius is well described by a compression model, Okuzumi et al's 
expression enables us to calculate the cross section, too. 

3. We further refined the compression model of SOS, by including the compression energy for 
the voids produced at the sticking of two aggregates. The refined model is also extended to 
non-equal-mass collisions in a reasonable way. With the refined model, we can accurately re- 
produce the evolution of both the gyration radius and the cross section of aggregates obtained 
by 808 from their early growth stage. The validity of the refined compression model for non- 
equal-mass collisions is also checked by additional numerical simulations of BCCA collisions. 
Although 808 considered only icy aggregates in the numerical simulation, our compression 
model would be also applicable to silicate aggregates by using a suitable value of E'roii- 

Our A^-body simulations of aggregate collisions and the refined compression model indicate that 
collisional compression is not so effective. As a result, dust aggregates (or initial planetesimal mate- 
rial) would have extremely low bulk densities, as suggested by 808. Okuzumi et al. (2012) showed 
that such extremely low bulk densities of aggregates accelerate their growth in protoplanetary disks 
and help their overcoming of the radial drift barrier against the planetesimal formation. However, 
solar-system bodies do not have such low densities at present. Dust aggregates (or planetesimals) 
should be compressed by other processes. A steady ram pressure due to the gas drag on aggregates 
and a self- gravity of sufficiently large aggregates would be candidates for aggregate compression, 
as indicated by 808 and Okuzumi et al. (2012). For relatively compact dust cakes (/) ~ O.lg/cm^) 
made of micron-sized silicate particles, compression is observed at a pressure > lOOPa (Blum and 
Schrapler 2004). However, for icy aggregates with very low bulk densities {p <C O.lg/cm^), com- 
pressive strength has not yet been measured. In future work, compression strength of very fluffy 
aggregates should be measured in numerical simulations and laboratory experiments. 

In the present study, we focus on the aggregates obtained at head-on collisions. At oblique 
collisions, the merged aggregates are elongated (WOT; Paszun and Dominik 2009). Although our 
model indicates inefficient compression at aggregate collisions, the effect of oblique collisions would 
further hinder compression. In future work, we should clarify the validity of our compression model 
in the case where oblique collisions are included. 



-15- 



The authors would like to thank Hiroshi Kimura, Tetsuo Yamamoto and Hiroshi Kobayashi 
for their valuable comments. We would also like to thank Takeshi Chigai for technical support with 
respect to the computer setup. This study was supported by a Grant-in- Aid from JSPS(22540242, 
22740299). 

REFERENCES 

Blum, J., & Schraplcr, R. 2004, Phys. Rev. Lett., 93, 115503 
Brauer, F., Dullemond, C. P., & Henning, Th. 2008a, A&A, 480, 859 
Brauer, F., Henning, Th., & Dullemond, C. P. 2008b, A&A, 487, LI 
Dominik, C, & Tielens, A. G. G. M. 1995, Phil. Mag. A, 72, 783 
Dominik, C., & Tielens, A. G. G. M. 1996, Phil. Mag. A, 73, 1279 
Dominik, C., & Tielens, A. G. G. M. 1997, ApJ, 480, 647 
Dullemond, C. P., & Dominik, C. 2005, A&A, 434, 971 
Goldreich, P., k Ward, W. R. 1973, ApJ, 183, 1051 
Goldstein, H. 1950, Classical Mechanics (Reading: Addison- Wesley) 

Giittler, C., Blum, J., Zsom, A., Ormel, C. W., & Dullemond, C. P. 2010, A&A, 513,A56 
Hayashi, C. 1981, Prog. Theor. Phys. Suppl., 70, 35 

Heim, L.-O., Blum, J., Preuss, M., & Butt, H.-J. 1999, Phys. Rev. Lett., 83, 3328 
Johnson, K. L., Kendall, K., & Roberts, A. D. 1971, Proc. R. Soc. London A, 324, 301 
Kempf, S., Pfalzner, S., k Henning, T. K. 1999, Icarus, 141, 388 
Krause, M., & Blum, J., 2004, Phys. Rev. Lett., 93, 021103 
Meakin, P. 1991, Rev. Geophys., 29, 317 

Min, M., Dominik, C., Hovenier, J. W., de Koter, A., Waters, L. B. F. M. 2006, A&A, 445, 1005 
Minato, T., Kohler, M., Kimura, H., Mann, I., & Yamamoto, T. 2006, A&A, 452, 701 
Mukai, T., Ishimoto, H., Kozasa, T., Blum, J., & Greenberg, J. M. 1992, A&A, 262, 315 
Nakagawa, Y., Nakazawa, K., & Hayashi, C. 1981, Icarus, 45, 517 
Ossenkopf, V. 1993, A&A, 280, 617 



- 16 - 



Okuzumi, S., Tanaka, H., & Sakagami, M. 2009, ApJ, 707, 1247 

Okuzumi, S., Tanaka, H., Kobayashi, H., k Wada, K. 2012, ApJ, in press ()arXiv:1204.5035p 
Ormel, C. W., Spaans, M., k Tielens, A. G. G. M. 2007, A&A, 461, 215 
Paszun, D., k Dominik, C. 2006, Icarus, 182, 274 
Paszun, D., & Dominik, C. 2008, A&A, 484, 859 

Safronov, V. S. 1969, Evolution of the Protoplanetary Cloud and Formation of the Earth and the 
Planets (NASA TTF-677; Moscow: Nauka) 

Stepinski, T. F., & Valageas, P. 1997, A&A, 319, 1007 

Sekiya, M. 1998, Icarus, 133, 298 

Suyama, T., Wada, K., k Tanaka, H. 2008, ApJ, 684, 1310 
Tanaka, H., Himeno, Y., k Ida, S. 2005, ApJ, 625, 414 

Wada, K., Tanaka, H., Suyama, T., Kimura, H., k Yamamoto, T. 2007, ApJ, 661, 320 

Wada, K., Tanaka, H., Suyama, T., Kimura, H., k Yamamoto, T. 2008, ApJ, 677, 1296 

Weidenschilling, S. J., k Cuzzi, J. N. 1993, in Protostars and Planets III, ed. E. H. Levy k J. I. 
Lunine (Tucson: University of Arizona Press), 1031 

Weidling, R., Guttler, C., Blum, J., k Brauer, F. 2009, ApJ, 696, 2036 

Youdin, A. N., k Goodman, J. 2002, ApJ, 620, 459 

Youdin, A. N., k Johansen, A. 2002, ApJ, 580, 494 

Youdin, A. N., & Shu, F. H. 2002, ApJ, 580, 494 

Zsom, A., Ormel, C. W., Dullemond, C. P., k Henning, Th. 2011, A&A, 534, A73 



This preprint was prepared with the AAS IATJhjX macros v5.2. 



-17- 



1mm 1cm 



CO 

CO „ 
< O) 

cr E 

LU ^ 
Q_ 

O ^ 

CO 

CO 
CO 

o 

DC 

o 




1 0"^ 1 
RADIUS [cm] 



Fig. 1. — Geometrical cross sections per masses for various types of aggregates. Constituent 
monomers are 0.1/im-seized icy spherical particles. The solid line shows the geometrical cross 
section of a compact sphere with the same density as constituent particles. The dotted line and 
the dashed line show the geometrical cross sections of BPCA and BCC A clusters, respectively. For 
BCCA and BPCA, we used empirical formulae bv lMinato et al.l (|2006l ). 
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Fig. 2. — Structure evolution of growing aggregates in the sequential collision simulations for various 
impact velocities Vimp and critical rolling displacements ^crit- Panel (a) shows the radius of gyration 
and panel (b) shows the density, respectively. The solid lines show the resultant aggregates in our 
simulations, and the dashed lines indicate BCCA clusters. Filled circles indicate the critical number 
of particles A'crit to start compression, as estimated from equation ([6]) with /3 = 0.5. 
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Fig. 4. — Geometrical cross sections of the resultant aggregates in the case of Vi^p = 0.54 m s~^ 
and = 2A as a function of the number of particles N. In panel (a), the dots show geometrical 
cross sections for 30 directions and the line shows the mean geometrical cross section. In Panel (b), 
thin gray lines show mean geometrical cross sections calculated for 30 resultant aggregates and the 
solid line shows the average value of them. We obtain evolution of the geometrical cross section 
with these two kinds of averaging. 
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Fig. 5. — Geometrical cross sections per unit mass of growing aggregates in the simulations for 
various impact velocities fimp and critical rolling displacements ^c- The solid lines show the ge- 
ometrical cross sections of the resultant aggregates in A^-body simulations, and the dotted lines 
indicate the geometrical cross sections of BCCA clusters. Collisional compression decreases their 
geometrical cross sections from the value of BCCA clusters. Filled circles indicate A'^crit given by 
equation dH) with /3 = 2.0. 
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Fig. 6. — The ratio of S to vrr^ of growing aggregates in the simulations for various impact velocities 
and critical rolling displacements. The solid lines show the resultant aggregates in our simulations, 
and the dotted line indicates that of BCCA clusters. We also plot the cross sections obtained from 
Okuzumi et al's expression (eq. [7]) with dashed lines. Okuzumi et al's expression reproduces the 
cross sections well even for compressed aggregates. 
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Fig. 7. — Same as Figure [6] but the prediction by the model of Paszun and Dominik (2009) is added. 
Only the case of = 8A is plotted. Okuzumi et al's expression reproduces the cross sections of 
compressed aggregates better than Paszun and Dominik's model. 
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Fig. 8. — Schematic illustration of three steps in the refined compression model. (a)Before a 
collision between two aggregates. (b)Step 1: At the moment of their stick. New voids are created 
in the merged aggregate. (c)Step 2: Compression of the new voids. (d)Step 3: Further compression 
of the merged aggregate. 
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Fig. 9. — Evolution of the gyration radius rg calculated with the refined compression model. The 
model curves are plotted with dashed lines. The parameters are set as 6 = 0.15 and b' = 36. The 
model curves agree well with the numerical results (solid lines) for all size range in all cases. 
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Fig. 10. — Evolution model of the cross section S obtained from equation ([7]) and the gyration radius 
calculated with the refined model (dashed lines). The refined compression model also reproduces 
the cross sections in the numerical simulations (solid lines) with the help of equation ([7]). 
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Fig. 11. — Same as Fig. [10] but the model cross sections are calculated with the refined model in 
the direct way, using the characteristic size rs- With the direct way, we can describe accurately 
the evolution of the cross section with the parameter b = 0.60. 
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Fig. 12. — Test by compression in non-equal-mass aggregate collisions. The results of N-hodj 
simulations of non-equal-mass BCCA collisions are plotted by circles and triangles for the mass 
ratios Mp/Mt = 1/16 and 1/4, respectively. Each data point is the averaged value obtained from 
10 runs of A^-body simulation for different types of initial BCCA clusters, as done by W08. The 
standard errors in the averaged values are also displayed. The predictions by the refined model are 
plotted by solid lines and the dashed line indicates the formula obtained by W08. 



